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On the July 2007 flare of the blazar 3C 454.3 
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ABSTRACT 

In July 2007, the blazar 3C 454.3 underwent a flare in the optical, reaching R ~ 13 on July 
19. Then the optical flux decreased by one magnitude, being R ~ 14 when the source was 
detected by the 7-ray satellite AGILE, that observed the source on July 24-30. At the same 
time, the Swift satellite performed a series of snapshots. We can construct the simultaneous 
spectral energy distribution using optical, UV, X-ray and 7-ray data. These shows that an 
increased 7-ray flux is accompanied by a weaker optical/X-ray flux with respect to the flare 
observed in the Spring 2005 by INTEGRAL and Swift. This confirms earlier suggestions about 
the behaviour of the jet of 3C 454.3. 
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1 INTRODUCTION 

Among the flat-spectrum radio quasar (FSRQ) class of blazars, 3C 
454.3 (z = 0.859) is one of the brightest sources. It has been de- 
tected significantly at almost all the available wavelength, from ra- 
dio to 7-rays (e.g., radio: Lister & Homan 2005, Homan & Lis- 
ter 2006; infrared: Bennett et al. 2003; optical: Villata et al. 2006, 
Fuhrmann et al. 2006, ultraviolet: Pian et al. 2005, X-rays: Tavec- 
chio et al. 2002, Pian et al. 2006, Giommi et al. 2006; 7-rays: Hart- 
man et al. 1999, Zhang et al. 2005, Nandikotkur et al. 2007). Its 
spectral energy distribution (SED) has the typical double-humped 
shape of blazars. In this case, the synchrotron peak is at infrared 
frequencies, while the inverse-Compton emission reach the max- 
imum at MeV-GeV energies (Fossati et al. 1998, Ghisellini et al. 
1998). 

In April-May 2005, 3C 454.3 literally exploded and had its 
historical maximum at optical wavelengths with i?=12.0 (Villata et 
al. 2006). This exceptional event triggered observations at X-ray 
energies with INTEGRAL (Pian et al. 2006) and Swift (Giommi et 
al. 2006). These detected the source up to 200 keV and in a bright- 
ened high energy state (Pian et al. 2006). This strong optical and 
hard X-ray activity was followed by a radio outburst with about 
one year delay (Villata et al. 2007). 

After a period of relative quiescence, the blazar underwent 
again an optical outburst in July 2007, although preliminary op- 
tical observation^ suggest a peak lower than that of 2005. The 



alert quickly spread all around the worlcj^ and triggered observa- 
tions with Swift and AGILE (Tavani et al. 2006). Although the lat- 
ter is still performing its verification phase, it was able to repoint 
towards the blazar on July 24-30, 2007 and to detect it (Bulgarelli 
et al. 2007). The flux measured by AGILE during the first part of 
monitoring, F = (3 ± 1) x 10~ 6 ph cm" 2 s _1 at E MOO MeV 
(Vercellone et al. 2007), is a bit less than 3 times the maximum 
[(1.25 ± 0.19) x 10~ 6 ph cm~ 2 s _1 ] and about 6 times greater 
than the average [(5.6 ± 0.4) x 10" 7 ph cm" 2 s" 1 ] flux observed 
by EGRET (Nandikotkur et al. 2007). 

Here we report about the analysis and assembling of the simul- 
taneous optical to 7-rays SED of 3C 454.3 in order to model this 
SED and test new ideas about how to obtain dramatic flux changes 
in specific bands without changing much the power carried by the 
jet, as detailed in the next section. 

We use H = 70 km s" 1 Mpc" 1 and « M = 0.3, f2 A = 0.7. 



2 SCIENTIFIC MOTIVATIONS 

During the huge multiwavelength flare observed in Spring 2005 
the optical flux increased by almost two orders of magnitude, while 
the X-ray flux increased by a factor ~10 with respect to the 2000 
level detected by BeppoSAX (Tavecchio et al. 2000). We can ask if 
this behavior corresponds to similar changes in the power carried 
by the jet or if, instead, the jet power is nearly constant, yet able 
to produce dramatic changes in specific energy bands. Pian et al. 
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(2006) could indeed describe the observed SEDs in 200(0 and 2005 
with minimal changes in the bolometric luminosity (and in the jet 
power), assuming that the dissipation region (where most of the 
radiation is produced) was inside the Broad Line Region (BLR) in 
2000, and outside it in Spring 2005. If the magnetic field is the 
same in the two epochs, despite the difference in the size of the 
emitting regions, then the power output is dominated by the Inverse 
Compton scattering of the broad line photons in 2000 (External 
Compton, EC), and by synchrotron and synchrotron Self Compton 
(SSC) emission in 2005. 

Then Katarzynski & Ghisellini (2007) proposed an alterna- 
tive model, in which the location of the dissipation region varies 
along the jet, and there is a link between this location, the com- 
pactness of the source, and the value of the bulk Lorentz factor F. 
The smaller this distance, the smaller T, the larger the magnetic 
field of the source. The physical bases of this model were i) the 
realization that the power lost in radiation must be a small frac- 
tion of the jet power, which has to carry most of it to the radio 
lobe regions (at least for powerful blazars); therefore ii) the kinetic 
and Poynting power of the jet are likely conserved, implying that 
the magnetic field B oc (Ti?) -1 , where R is the (transverse) size 
of the jet; iii) in the internal shock model for blazars the dissipa- 
tion is due to the collision of two shells moving with different F: 
the collision occurs at a distance which is a factor F 2 larger than 
the initial shell separation. This introduces a link between com- 
pactness of the source and F; iv) even if the dissipation is not due 
to internal shocks, we may have a similar relation between V and 
compactness if the jet is still accelerating when dissipation occurs. 
Following these ideas, Katarzynski & Ghisellini (2007) could re- 
produce the 2000 and 2005 states of 3C 454.3 with a jet of exactly 
the same power, by assuming different F (11 in 2000 and 6.25 in 
Spring 2005). In this model the observed bolometric luminosity is 
nearly constant (except for extreme changes in F, resulting in ex- 
treme changes of the Doppler factor and thus in the Doppler boost- 
ing of the observed radiation). When dissipation occurs with large 
r the seed external photons are seen largely boosted, while the B- 
field is relatively weaker, resulting in a EC dominated spectrum, 
where most of the power is emitted the MeV-GeV band. On the 
other hand, for smaller F, the 73-field is larger, the external ra- 
diation is seen less boosted, resulting in a synchrotron and SSC 
dominated spectrum, with much less contribution of the EC. 

It is therefore clear that the main diagnostics of this model is 
the level of the MeV-GeV flux during flares: if both the optical 
(synchrotron) emission and the MeV-GeV (EC) emission increase, 
then this flags a real increase of the jet power, but if the increase of 
the MeV-GeV flux corresponds to the same (or even to a fainter) 
optical flux, then this flags an emission episode by a jet of nearly 
constant power but larger F. Vice-versa, if the MeV-GeV radiation 
decreases during an optical flare, this flags a jet of constant power 
dissipating closer to the black hole and with a smaller F. 

We hope that these simple trends will be verified in the era of 
simultaneous GLAST and Swift observations, able to look for the 
crucial part of the SED (optical to GeV bands). But anticipations of 
this, for the brightest blazars, are possible using AGILE and Swift 
observations. This is exactly what happened for 3C 454.3 in July 
2007. 



3 The 2000 SED actually includes also data taken at previous epochs. We 
will hereafter refer to it as "2000 SED" because the modelling of the SED 
is critically driven by the BeppoSAX spectrum, acquired in 2000. 



3 DATA ANALYSIS OF THE SWIFT OBSERVATIONS 

Although Swift observed 3C 454.3 several times in July and August 
2007, we analysed data from the only observation simultaneous 
with AGILE, i.e. that performed on July 26. The data from the three 
instruments onboard Swift (Gehrels et al. 2004) have been pro- 
cessed and analysed with HEASoft v. 6.3.1 with the CALDB 
release of July 30, 2007. 

Data from the X-ray Telescope XRT (0.3 — 10 keV, Burrows 
et al. 2005) were analysed using the xrtpipeline task, set for 
the photon counting mode and having selected single to quadruple 
pixel events (grades 0—12). The output data were rebinned in order 
to have at least 30 counts per energy bin. The final spectrum is best 
fitted (x 2 = 45.9 for 43 de grees of freedom) with a power law 
with absorption in excess of the Galactic value (N^^ 1 — 7.24 x 
10 20 cnT 2 , Kalberla et al. 2005). Formally, this additional at 
z = 0.859 hasavalue of 1.9li'J X 10 21 cm -2 . The power law has a 
photon index of 1.47±0.11 and a flux of 6.1 x 10 -11 erg cm -2 s _1 
in the 2 — 10 keV energy band. 

The additional iV^ may be not physical, but it is used to fit 
satisfactorily the observed deficit of photons at low energies. A fit 
with a broken power-law model gives multiple minima for the low 
energy photon index, making it impossible to constrain efficiently 
the parameters. The photon deficit at low energies has a gradual 
decrease, resembling a low-energy cutoff, and the additional ab- 
sorption is the best statistical way to mimic it. Therefore, when in 
the following we refer to the unabsorbed X-ray spectrum, we mean 
the spectrum corrected only for the Galactic absorption. 

Data from the UV telescope UVOT (Roming et al. 2005) were 
analysed by using the uvotmaghi st task with a source region of 
10", while the background was extracted from a source-free circu- 
lar region with radius equal to 60" (it was not possible to use an 
annular region, because of a nearby source). We added a 10% error 
in flux (corresponding to about 0.1 magnitudes) to take into ac- 
count systematic effects. The observed magnitudes are UVW1 = 
14.7 ± 0.1, UVM2 = 14.9 ± 0.1, and UVW2 = 15.1 ± 0.1. 

Instead, for the hard X-ray detector BAT (optimised for the 
15 — 150 keV energy band, Barthelmy et al. 2005), it is necessary 
to perform long exposure observations in order to get useful signal 
from extragalactic sources. Therefore, data from all the July obser- 
vations (performed on days 11, 12, 13, 14, 17, 20, 26, 28, 30) were 
binned, cleaned from hot pixels and background, and deconvolved. 
The intensity images were then integrated by using the variance as 
weighting factor and the total exposure was 15 ks. The flux upper 
limits for a 3a detection - already corrected for systematics - mea- 
sured at the position of 3C 454.3 are 2.1 x 10" 10 erg cm -2 s" 1 (28 
mCrab) and 3.6 x 10~ 10 erg cm" 2 s" 1 (38 mCrab) in the 20-40 
keV and 40-100 keV energy bands, respectively. 



4 MODELLING THE SED 

We adopt a simple, one-zone, homogeneous synchrotron and in- 
verse Compton model to the different states of 3C 454.3. The model 
is fully described in Ghisellini, Celotti & Costamante (2002), and 
we briefly outline here the general assumptions. 

The source is a cylinder of cross sectional radius ii = ipz, 
where z is the distance from the apex of the jet, assumed to be a 
cone of semiaperture angle ip ~ 0.1. The width AR', as measured 
in the comoving frame, is assumed to be within the range R/F ^ 
AR' ^ R. The magnetic field B is homogeneous and tangled; the 
blob moves with a bulk Lorentz factor F and the viewing angle 
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Figure 1. SEDs of 3C 454.3 at different epochs. Top panel: the SED in 
2000, corresponding to the BeppoSAK observations discussed in Tavec- 
chio et al. (2002). The other points are not simultaneous, (see Pian et al. 
2006 and references therein). Mid panel: the SED during the huge optical 
flare in 2005, as described in Pian et al (2006) and Giommi et al. (2006). 
Different colors correspond to different sets of strictly simultaneous data. 
Bottom panel: the SED on July 26, 2007, as observed by AGILE, and Swift. 
We have also added the optical flux in the R band, taking the average be- 
tween the magnitude of July 24 and July 29 measure by the Tuorla obser- 
vatory. The optical and X-ray fluxes are corrected for Galactic extinction 
(Ay = 0.355). The solid and dashed lines correspond to our modelling 
(parameters in Tab. [TJ. The dotted line is the contribution from the accre- 
tion disk (assumed to be a simple black-body). 

is 9. The external radiation is produced at a fixed radius, that can 
be identified with the radius .Rblr of the broad line region (BLR; 
Sikora, Begelman & Rees 1994), reprocessing ~10% of the disk 
luminosity. This is likely an approximation of the real amount of 
external radiation, since, besides the contribution of the BLR, other 
processes can contribute (scattering by ionised intercloud plasma, 
synchrotron radiation "mirrored" by the clouds and/or the walls of 
the jet, Ghisellini & Madau 1996, and also direct radiation from the 
accretion disk, Dermer & Schlickeiser 1993; Celotti, Ghisellini & 
Fabian, 2007). 

The particle distribution is the result of injection (for a finite 
time, taken to be equal to AR'/c in the comoving frame) and cool- 
ing. We calculate the random Lorentz factor 7 CO oi at which the par- 



ticles cool in one light crossing time. If the particles are injected 
between -ytnj and 7max with a power law distribution of slope s, in 
the "fast cooling" regime (7 coo i < 7mj), we have N(j) cx 7~ p 
between 7i n j and 7 ma x, (where p = s + 1 is the injection slope 
increased by one unit) and AT (7) oc j~ 2 between 7 CO oi and 7^. In 
the models presented here, this is always the case, since 7 coo i is of 
order unity in all cases. 

Based on the above assumptions, the SED modelling yields 
the results shown in Fig.Q] where we report model curves for four 
states of the source. The upper panel shows the SED correspond- 
ing to the BeppoSAX observations in 2000, with which we asso- 
ciate the (not simultaneous) COMPTEL, EGRET and optical data 
(Tavecchio et al. 2002; Pian et al. 2006; Hartman et al. 1999; Zhang 
et al. 2005). The mid panel shows two states in 2005, presented in 
Pian et al. (2006) and Giommi et al. (2006), for which we have 
simultaneous optical, soft and hard X-ray INTEGRAL data. The 
bottom panel shows the state in July 26, 2007, as observed by Swift 
and AGILE, together with an optical point (in the R band) taken 
during the multiwavelength campaign of this source by the Tuorla 
observatory. 

The input parameters used for these models are listed in the 
upper part of Tab.Q] In the lower part we report the derived Doppler 
factor S — l/[r(l — /3cos9)] and the power carried by the jet 
in the form of magnetic field (Lb), cold protons (L p ), relativistic 
electrons (L c ), and produced radiation (L ra d), defined as: 

L B = 7yR 2 T 2 cU b 

L p = -kR 2 Y 2 c J N (-f)m p c 2 d'y 

L = ttR 2 V 2 c J N (■y)'ym c c 2 d'y 

Lrad = 7vR 2 F 2 cULd (1) 

where Ub and (7 r ' ad are the magnetic and the radiation energy den- 
sity measured in the comoving frame, respectively, and where we 
have assumed one proton per relativistic electron. 

The choice of the parameters is not unique in SSC+EC mod- 
els, especially when the high energy (7-ray) peak is not observed, 
as in 2005. The adopted choice reflects the following considera- 
tions: i) we have used the low optical state in 2000 (top panel) to 
fix the contribution of the accretion disk luminosity; ii) we have 
tried to minimise the total power budget; iii) we have tried to de- 
scribe the different states of the source with a minimal change of 
the power carried by the jet. 

We find that the July 2007 flare and the state in 2000 can 
be reproduced by very similar parameters (same radius and bulk 
Lorentz factors, slightly different magnetic field and shape of the 
particle distribution), the main difference being the amount of in- 
jected power (-t/Jnj going from 7 x 10 43 to 3 x 10 44 erg s _1 ). For 
the two states in 2005, instead, the SED can be reproduced by a 
more compact source, with a smaller T, corresponding to a smaller 
Doppler factor 5. The magnetic field is greater (B — 15-35 G), 
and the injected power is slightly greater than for the 2007 state. 

The model can reproduce the huge variations observed espe- 
cially in the optical (between 2000 and 2005 there is a difference of 
a factor ~ 100) and in the 7-ray band (the 2007 state is one order of 
magnitude larger than in 2000). This is achieved with much smaller 
variations of the power transported by the jet, as discussed below. 

Note that in all cases the emission drops above 50-100 GeV. 
This is due to the robust limit to the maximum emitted 7-ray en- 
ergy provided by the decrease of the Klein-Nishina cross section. 
The maximum energy can be estimated to be /ii^nax — 6S/vt,i5 
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Table 1. For all models we have assumed a viewing angle 9 = 3.5° and a 
bolometric luminosity of the accretion disk Ldisk = 3 X 10 46 erg s _1 . We 
assume that the 10% of the disk luminosity is reprocessed by the BLR and 
re-emitted as broad lines at a distance 7?blr = 3 X 10 17 cm for all models. 
The minimum variability timescale t V ar is defined as i va r = R/ (ci5). 
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Figure 2. X-ray data of July 2007 with, superimposed, the result of our 
model assuming that the minimum energy of the emitting electrons corre- 
sponds to 7 m j n = 1,2 and 3, as labelled. Only for 7 m i n = 1 we obtain a 
reasonable fit. 

GeV, where Ut,i5 is the frequency of the target soft photons (as 
seen in the observer frame) in units of 10 15 Hz. Internal (i.e. inside 
the source) photon-photon absorption, instead, plays no role in the 
models presented here. 



5 DISCUSSION 

Our modelling of the SED of 3C 454.3 in 2000, 2005 and 2007 is 
in qualitative agreement with the expectations of the jet radiation 
model proposed by Katarzynski & Ghisellini (2007), where the 
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Fi gure 3. The location of the different states of 3C 454.3 in the 7peak — 
(Ur + Us) plane. Adapted from Foschini et al. (2006). 

relative importance of the synchrotron and SSC luminosity (with 
respect to the EC one) is controlled by the value of T, which is as- 
sociated with the compactness of the source. The relation between 
r and source compactness relies either on the assumption that the 
jet is still accelerating when it dissipates, or on the assumption that 
the dissipation is the result of internal shocks between two shells 
(Spada et al. 2001), occurring at a distance ~ T 2 times the initial 
separation: in both cases we have that emitting regions with smaller 
r are closer to the black hole, and thus more compact. The mag- 
netic field would be stronger, implying larger synchrotron and SSC 
luminosities with respect to the EC one. It is then possible to have 
dramatic variations in specific frequency bands even for jets carry- 
ing the same amount of power. 

Indeed, the power in bulk motion of the relativistic electrons 
(see Tab. [TJ varies no more than a factor 6.3 for the 4 different 
states we have considered, while the Poynting flux varies merely 
by a factor 3. The largest change of the jet power is due to the pro- 
ton component (factor 8.6), assuming that there is, for all states, 
one proton per relativistic electron, and that the electron distribu- 
tion extends down to 7 m i n = 1. The latter assumption is usually 
very difficult to test, since low energy electrons emit synchrotron 
radiation in the self-absorbed regime (which is both insensitive to 
the low-energy cutoff and practically impossible to observe) and 
SSC radiation at relatively low frequencies, where the flux is dom- 
inated by the synchrotron emission of electrons of higher energies. 
The best opportunity to see the radiation emitted by electrons of 
low 7 is thus when the SSC contribution is not contributing much 
to the emission in the X-ray band, which is instead dominated by 
EC radiation. In fact in this band the EC is produced by low energy 
electrons. This is the case for the 2000 and 2007 states of 3C 454.3. 
Fig. |2] shows the predicted spectrum varying 7 m in from 1 to 3. It 
shows that even tiny variations have a large effect on the spectrum, 
which becomes much harder than observed when 7 m i n > 1. We 
then conclude that the electron distribution must extend down to 

7min ~ 1. 

Consider now the value of L ra( j in comparison with the other 
powers. It is of the order of, but in some case larger, than L e and 
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Lb- This shows that the power spent by the jet to produce the ra- 
diation we see cannot be provided by the magnetic field or by rel- 
ativistic electrons only. There must be another form of power, able 
to re-furnish the electrons of the energy lost for producing the radi- 
ation. The simplest possibility is that the jet is baryon loaded, and 
that L p dominates the jet power. In our case (assuming one proton 
per electron) we can estimate the radiative efficiency (i.e. L ra d /L p ) 
which is between 4 and 12 per cent. This is in excellent agreement 
with the idea that the jet must be radiatively inefficient, since it has 
to carry most of its power to the radio lobes (see also Celotti & 
Ghisellini 2007). 

Finally, we can see where the different states of 3C 454.3 are 
located in the plane 7 pca k-f > where 7 pea k is the random Lorentz 
factor of the electrons emitting at the two peaks of the SED, and 
U — Ub + £/rad is the radiation+magnetic energy density as seen 
in the comoving frame. 

In this plane blazars form a well defined sequence (see Ghis- 
ellini et al. 1998; 2002): low power BL Lacs follow a relation 
7 P cak oc (7 _1 and more powerful blazars (FSRQ with broad emis- 
sion lines) tend to lie on a 7 pca k oc U~ ' 6 relation. As Fig.[3]shows, 
the 4 states of 3C 454.3 considered here are somewhat scattered, but 
they are within the total scatter of the relation in this high U, small 
7 pca k part. Note also that in this plane the 2000 and 2007 states (of 
strong 7-ray flux) are identical. 

We can conclude that the July 2007 flare of 3C 454.3 gives 
support to the idea that dissipation site in relativistic jets changes, 
and that it is likely that, when it is located closer to the black hole, 
the emitting region is more compact, has a smaller bulk Lorentz 
factor, and a greater magnetic field, although the power in Poynting 
flux needs not to change, since a larger B-field is compensated by 
a smaller size and bulk Lorentz factor, making Lb nearly constant. 
Observationally, these characteristics translate in dramatic varia- 
tions of the synchrotron peak flux and of the EC high energy com- 
ponent, and in an inverse correlation between the synchrotron and 
EC components. The more compact the source is, the more intense 
its synchrotron flux, at the expenses of the EC emission. The lat- 
ter becomes fainter for the weaker Doppler boosting of the exter- 
nal seed photons in the comoving frame, due to the smaller bulk 
Lorentz factor. The flux variations at high energies are however mit- 
igated by the SSC emission, which is enhanced for more compact 
emission regions. This component can dominate in the X-ray band, 
and can be recognized by a softer and more intense X-ray spectrum 
with respect to the spectrum due to the EC process. Since in this 
case the production of high energy 7-rays is reduced, the pair pro- 
duction process, due to photon-photon interactions, is reduced also, 
and we can evade the problem of the reprocessing of high energy 
7-ray luminosity into the X-ray band, which inhibits the produc- 
tion of GeV radiation in too compact sources. If different states are 
indeed related to different compactness of the source, one might 
also expect different variability timescales. However (and unfortu- 
nately) this is mitigated by the different Doppler factors, which are 
smaller for more compact emitting regions, resulting (see Table 1) 
in minimum variability timescales f var = R/(cS) that are almost 
equal. 
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